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Abstract 

We use a density matrix renormalization group method to study quantita- 
tively the phase diagram of a one-dimensional extended Hubbard model at 
half- filling by investigating the correlation functions and structure factors. We 
confirm the existence of a novel narrow region with long-rang bond-order- wave 
order which is highly controversial recently between the charge-density-wave 
phase and Mott insulator phase. We determined accurately the position of the 
tricritical point Ut — 7.2t, Vt — 3.746t which is quite different from previous 
studies. 
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Much effort has been devoted to understanding the effects of strong electronic correla- 
tions on quasi-one-dimensional systems. These materials including conjugated polymers, Cu 
oxides and Ni halides exhibit rich phase diagrams and display huge ultrafast optical non- 
linearity which point to promisinig optoelectronic applications [1-4]. The one-dimensional 
(ID) extended Hubbard model (EHM) with the nearest neighbour repulsion V in addition 
to the on-site repulsion U is a standard minimal model that can describe these rich phys- 
ical properties. Although this model has been investigated for more than two decades [5], 
its ground-state phase diagram is still controversial [6-16]. In the limit V = [17], the 
ground state is in the Mott insulating state where the spin sector shows quasi-long-range 
order of spin density wave. With use of weak-coupling renormalization group analysis [5], a 
continuous phase transition is obtained at U = 2V from charge-density-wave (CDW) to spin- 
density-wave (SDW) phases. Strong-coupling calculations using second order pertubation 
theory [5] gave a first-order phase transition also at U — 2V. This means that there exists 
a crossover between these two phase transitions at a tricritical point in the intermediate- 
coupling regime. This picture was supported by both numerical [6,7,10] and analytical [7-9] 
studies and had been regarded for a long time as the complete phase diagram of the ID 
EHM at half filling. However, the precise location of the tricritical point where the nature 
of the transition changes has been a subject of much investigations and remains uncertain 
[6-16]. 

Recently, by using level crossings in the excitation spectra obtained by exact diagonal- 
ization, Nakamura [11] pointed out that a novel spontaneously dimerized phase which is also 
called bond-order-wave (BOW) phase exists for small to intermediate values of U and V in 
a narrow strip between CDW and SDW phases. This remakable proposal was confirmed by 
stochastic series-expansion Monte Carlo calculations [13,15] and reformulated weak-coupling 
field theory [12] where higher-order terms were included but was questioned by a high-level 
density-matrix renomalization group (DMRG) investigation [14] which shows the existence 
of a BOW phase only on a segment of a first-order SDW-CDW phase boundary. On the 
other hand, Zhang [16] cast doubts on the accuracy of the DMRG calculation and predict 
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that a SDW-CDW phase boundary smoothly connects the weak- and strong-coupling limits. 

In order to clarify these apparent contradictions, we perform an extensive numerical study 
of ID EHM at half-filling in the vicinity of U = 2V using the DMRG [18] technique. In 
contrast to the previous DMRG studies [10,14,16], we mainly calculate the BOW correlation 
functions and structure factors which are the most direct evidences for the long-range BOW 
phase. Carefully using various finite-size scaling skills for different physical quantities defined 
below, we confirm that the recently discovered BOW phase [11] of the ID EHM does have a 
finite extent in the (U, V) plane, but its width is much smaller as has been predicted by the 
Monte Carlo method [13]. The position of the tricritical point where the BOW phase vanishes 
is determined accurately. The value is much higher than any previous results [6-8,10,12-14] 
while suprisingly it is consistent with the point given in ref. [14] where Jeckelmann argued 
that the BOW phase exists up to an upper limit which is smaller than 8t. However the point 
is not identical to the tricritical point. This can be understood as the result of increasing 
frustration in the spin degrees of freedom [9,14]. Furthermore we find that the first-order 
phase transition takes place for U < 2V rather than U > 2V obtained from reformulated 
weak-coupling field theory [12]. Using the same methods described below for U up to lQt 
with the interval of AU = 0.2t, we quantitatively determine the phase diagram of the ID 
EHM at half-filling for the first time. 

The Hamiltonian of the ID EHM is defined as follows 

H = -tJ2 Bi,i+i + n Wl + V Y1 n ' n '+ 1 ' w 
i i i 

where = c^ a c\^ — |, ni = + nil, c ta denotes the creation operator of an electron at 

the Ith site with spin a and the bond-charge density operator B^i +1 is 

Bl,l+i = i C ta C l+l,a + cj_ 1)ff ci )<r ) • (2) 

cr 

As is well-known, DMRG is a very accurate numerical method for the ground-state 
properties of a ID quantum system with short-range interactions [18]. Here we applied the 
finite-size DMRG algorithm with open boundary conditions to study the Hamiltonian (1) 
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at half-filling. This method allows us to probe directly correlation functions and structure 
factors associated with SDW, CDW and BOW in the ground state. Lattices up to 512 sites 
were used in our studies. The largest number of states kept in the calculation was m = 512 
per block. The hopping integral t is set to 1 as the energy unit. The weight of the discarded 
states was typically about 1CT 7 — 1CT 14 depending on whether the system is in its critical 
state or not in the final sweep. The convergence tests as functions of number of states kept 
and system size were carefully performed. We checked our DMRG calculations against exact 
numerical results for noninteracting (U — V — 0) chains (up to 512 sites) and results from 
exact diagonalization for interacting (U ^ 0, K / 0) chains (up to 14 sites). Excellent 
agreement was found in both cases. Furthermore, for a given chain length, we extrapolate 
our DMRG results to the limit of vanishing discarded weight. When interactions are turned 
on, there exist finite excitation gaps on finite chains, so the accuracies of all quantities we 
calculated are not worse than those of the noninteracting case. Thus, numerical errors in 
our work could be safely estimated to be smaller than 10~ 4 . 

First, we will introduce a method which can be simply used to determine the phase 
boundaries of BOW phase directly. The BOW structure factor is given by 

Sbow{q) = — e tqr ((Bij + iBi +r j +r+ i) 

lr 

- (Bi,i+i) (B l+r j +r+1 )) . (3) 

As pointed out by Nakamura [11], phase transitions on the CDW-BOW and BOW-MI phase 
boundary are of Gaussian type and Kosterlitz-Thouless (KT) type respectively, the real 
space staggered bond fluctuation correlation function falls off algebraically as 

(-l) r «£ J)J+1 £ J+r)J+r+1 ) - (B hl+1 ) (B l+r , l+r+1 )) ~ r"" (4) 

Away from phase boundaries, this quantity falls off exponentially. Thus for a finite-size 
system, the Sbow ij) is expected to reach a maximum at the critical points. Fig. 1 shows 
the results of the Sbow (tt) for different system sizes with U/t = 4.0 and 1.5 < V/t < 2.5. 
As expected, the Sbow i 71 ) peaks twice for all the different system sizes we calculated, 
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clearly indicating that there exist two phase transitions. The inset of Fig. 1 shows a linear 
extrapolation of these two critical values with the inverse of the chain length 1/L. We find 
that the larger the system is, the smaller the BOW phase becomes. Nevertheless, the BOW 
phase remains finite at the thermodynamic limit. Compared with other recent numerical 
methods [11,13,14], we find this to be the most accurate way to locate the BOW phase 
boundaries. 

Then we will further prove the existence of a finite-width BOW phase established by our 
method. The most direct evidence is the BOW correlation function 

Cbow (r) = (-l) r {\ (Bi, i+1 B l+r , l+r+1 ) - B 2 ) (5) 

i 

where B = jr £^ I n Fig. 2 (A), we show the staggered BOW correlation functions 

with increasing nearest neighbour repulsion V at U/t = 4.0. To avoid boundary effects, 
we only perform the average in (5) over 256 sites in the middle of the 512-site system. 
The results indicate that there exist three different phases in model (1) since the staggered 
BOW correlation functions show three distinct type of behavior as r increases: (i) it decays 
exponentially at V/t = 1.5, indicating that the system has no BOW order; (ii) it converges to 
a nonzero constant, indicating that the system is in the long-range BOW phase at V/t = 2.1; 
(iii) it decays as 1/r at V/t = 2.5, indicating that the system is in the Mott insulating 
phase. Besides the data we show in Fig. 2 (A), we also compare our DMRG results with 
Quantum Monte Carlo calculation [13] at U/t = 4.0, V/t = 2.14 for L = 128 and 256. The 
amplitudes of the oscillations at large distances are the same. The BOW order parameter 
in the thermodynamic limit 

A BOW = hm Ij](-1)'(£ M+1 > (6) 

I 

can be obtained by fitting ttibow (L) (=y^2 CBow( r )/L) with a second-order polynomial 
in 1/L since m B ow (L) — > A B ow for L — > oo. Fig. 2 (B) shows such extrapolations. We find 
that m 2 BOW (L) approaches zero when V/t = 1.5 and 2.5 but remains finite when V/t = 2.1. 
The strength of dimerization is consistent with Quantum Monte Carlo calculation [15] while 
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another recent DMRG investigation [14] has not detected this rather strong BOW order. 

In order to give more convincing evidence, we did another finite-size analysis in the 
vicinity of these two phase transitions. Let us start from the first phase transition at U — U c . 
Fig. 2 (C) presents plots of \b.[Sbow (tt)] versus ln[L] lorU/t — 4.0 and three different values 
of V/t around the first critical point. Data points for V/t = 2.16 indeed fall on a straight 
line, indicating critical scaling for the BOW fluctuations. At the other two points V/t — 2.15 
and 2.17, data points behave nonlinearly due to the exponential decay term. Applying the 
same finite-size analysis to the CDW structure factor 

Scdw (q) = j^2 elQr (( n i n i+r) - (ni) (n l+r )) , (7) 

Ir 

we can also explore the nature of the first phase transition. The linear behavior of 
ln[ScDW (tt)] around V/t = 2.16, shown in Fig. 2 (D), confirms the vanishing of the charge 
gap at the first phase transition point. 

Next we determine the location of the phase transition of KT type. As argued by 
Nakamura [11], it is a quantum phase transition of the KT type. This makes it difficult to 
determine the phase boundary directly from the behavior of Ssdw ( n ) 

SsDw(q) = jJ2 et9r ( S l S ^r) (8) 
Ir 

due to the finite-size effects. Instead, we apply an indirect method, used by Sengupta et al. 
[13,19], to confirm the second phase transition. It is well known [8] that if the ground state 
of a ID system is spin-gapless, the spin-spin correlation falls algebraically with exponent 
equal to 1. It has been further shown [20] that in the spin-gapless phase Ssdw (^) — ^ I/ 71 " 
as q — > whereas in the spin-gapped phase Ssdw (<?) /<? 0. Even a very small spin gap 
can be detected in this way, since it is in practice sufficient to see the 7iSsdw (<?) /q decay 
below 1 for small q to conclude that a spin gap must be present. From Fig. 2 (E) we can 
see the behavior of ttSsdw (<?) /q for U/t = 4.0 and different values of V/t. In the gapless 
region, logarithmic corrections [21] make it difficult to observe the approach to 1 as q — > 0. 
In analogy with spin systems [22], we expect the leading logarithmic corrections to vanish 

6 



at the point where spin gap opens and therefore exactly at the critical point there should 
be a clear scaling to I. Based on results shown in Fig. 2 (E), we estimate the BOW-MI 
boundary to be at V/t = 2.0 ± 0.01 at U/t = 4.0 which is consistent with the results shown 
in Fig. 1. The Fig.2 (F) further provides evidence for the transition from the spin-gapped 
state, as identified by the exponential decay of the staggered SDW correlation function, to 
the spin-gapless state, as characterized by the 1/r-decay of the SDW correlation function 



Now Let us determine the tricritical point. As we mentioned above, Nakamura [11] first 
argued that there exists an extended BOW phase for couplings weaker than a tricritical 
point, which is supported by numerical [13] and analytical [12] studies where the tricritical 
point is given as U t ~ 4.76t and U t ~ 5t respectively. However, recent DMRG calculation 
[14] disagreed this picture and pointed out that the BOW phase exists only on the first-order 
SDW-CDW phase-transition line for intermediate coupling U starting from the tricritical 
point U t ~ 3.7t to an upper limit which is smaller than 81 In order to give a more convincing 
result, we have studied the phase transitions in a wide range of coupling U starting from 
2t to 9t with the interval of AU = 0.2t as carefully as we displayed above. We find that 
the width of the BOW phase will shrink to at U ~ 7.2t in the thermodynamic limit 
and the extrapolation of the BOW order parameter will go to for all the couplings V at 
U > 7.2t which indicate the vanishing of long-range BOW phase and the merging of these 
two continuous phase transitions into one. Furthermore, we also calculate the V dependence 
of the CDW order parameters m C DW (L) (the definition is similar to m BO w (L)) across the 
phase boundary for U up to 16. From the inset of Fig. 3 we can see that the characteristics 
of a second-order transition for U = 6t and a first-order transition for U — 9t are indeed 
quite apparent. The points are obtained from calculating the CDW order parameters for 
the chain length L = 192, 256, 320, 512 and then extrapolating to the thermodynamic limit 
L — > oo. The transition point V c ~ 4.143 for U — 8t and V c ~ 6.116 for U = 12t are 
consistent with previous results [13,14]. 




(9) 
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The phase diagram obtained from DMRG calculations described above is shown in Fig. 
3. The BOW phase exist from weak coupling up to tricritical point where the two continuous 
transition line meet. With the increasing of U and V along the line U = 2V, the BOW 
phase first expands and then shrinks up to the tricritical point. Beyond this point the 
BOW phase disappears and the first order transition from the CDW to SDW phase can 
be obtained. This phase diagram is similar to the schematic one in ref. [13] except that 
the width of BOW phase become much smaller and extend to much larger coupling. The 
first-order phase transition takes place for U < 2V rather than U > 2V reported in ref. [12]. 
Our result does not agree with the location and the width of the BOW phase given in ref. 
[14] while the upper limit of the BOW phase is consistent with each other. 

In conclusion, we have studied ID half-filled EHM using the DMRG method. The quanti- 
tative phase diagram is obtained by investigating correlation functions and structure factors 
for the first time. The novel spontaneously dimerized phase does have finite width in the 
(U, V) plane. Our estimate for the tricritical point where the BOW phase vanishes and the 
nature of the transition changes is U t /t ps 7.2, V t /t ps 3.746. The extension of dimerization 
to rather strong coupling can be understood as the results of increasing frustration in the 
spin degrees of freedom [9,14]. 

It is our pleasure to acknowledge helpful discussions with Dr. P. Thalmeier. 



8 



REFERENCES 

[1] Conjugated Conducting Polymers, edited by H. Kiess ( Springer- Verlag, Berlin, 1992). 

[2] T. Ogasawara, M. Ashida, N. Motoyama, H. Eisaki, S. Uchida, Y. Tokura, H. Ghosh, 
A. Shukla, S. Mazumdar, and M. Kuwata-Gonokami, Phys. Rev. Lett. 85, 2204 (2000). 

[3] H. Kishida, H. Matsuzaki, H. Okamoto, T. Manabe, M. Yamashita, Y. Taguchi, and Y. 
Tokura, Nature (London) 405, 929 (2000). 

[4] S. Iwai, M. Ono, A. Maeda, H. Matsuzaki, H. Kishida, H. Okamoto, and Y. Tokura, 
Phys. Rev. Lett. 91, 057401 (2003). 

[5] V. J. Emery, in Highly Conducting One- Dimensional Solids, edited by J. T. Devreese, 
R. Evrard, and V. van Doren (Plenum, New York, 1979). 

[6] J. E. Hirsch, Phys. Rev. Lett. 53, 2327 (1984). 

[7] J. W. Cannon and E. Fradkin, Phys. Rev. B 41, 9435 (1990); J. W. Cannon, R. T. 
Scalettar, and E. Fradkin, ibid. 44, 5995 (1991). 

[8] J. Voit, Phys. Rev. B 45, 4027 (1992). 

[9] P. G. J. van Dongen, Phys. Rev. B 49,7904 (1994). 

[10] G. P. Zhang, Phys. Rev. B 56, 9189 (1997). 

[11] M. Nakamura, J. Phys. Soc. Jpn. 68, 3123 (1999); Phys. Rev. B 61,16377 (2000). 

[12] M. Tsuchiizu and A. Furusaki, Phys. Rev. Lett. 88, 056402 (2002); cond-mat/0308157. 

[13] P. Sengupta, A. W. Sandvik, and D. K. Campbell, Phys. Rev. B 65,155113 (2002). 

[14] E. Jeckelmann, Phys. Rev. Lett. 89, 236401 (2002); 91, 089702 (2003). 

[15] A. W. Sandvik, P. Sengupta, and D. K. Campbell, Phys. Rev. Lett. 91, 089701 (2003). 

[16] G. P. Zhang, Phys. Rev. B 68, 153101 (2003). 

9 



[17] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20,1445 (1968). 

[18] S.R. White, Phys. Rev. B 48, 10345 (1993). 

[19] P. Sengupta, A. W. Sandvik and D. K. Campbell, Phys. Rev. B 67, 245103 (2003). 

[20] R. T. Clay, A. W. Sandvik and D. K. Campbell, Phys. Rev. B 59, 4665 (1999). 

[21] S. Eggert, I. Affleck and M. Takahashi, Phys. Rev. lett. 73, 332 (1994). 

[22] S. Eggert, Phys. Rev. B 54, 9612 (1996). 



10 



FIGURES 

FIG. 1. Behavior of Sbow{^) across the phase boundary for U/t = 4.0. The inset shows a 
linear extrapolation of the critical values V c and V s with the inverse of chain length 1/L. 

FIG. 2. Finite-size analyses for BOW, CDW, and SDW correlation functions and structure 
factors in the vicinity of the two phase transitions at V c and V c for U/t = 4.0. 

FIG. 3. Phase diagram of the half-filled extended Hubbard model. The tricritical point is at 
(U t ,V t ) ~ (7.2i, 3.746f). The inset shows V dependence of the CDW order parameters mcDW 
across the phase boundary for U = 6t and U = 9t in the thermodynamic limit. 
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